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AN ANALYSIS OF BUZZING IN SUPERSONIC RAM JETS BY A 


MODIFIED ONE-DIMENSIONAL NONSTATIONARY 


WAVE THEORY 
By Robert L. Trimpi 


SUMMARY 


Experimental instantaneous pressure records of simulated ram- jet 
models without heat addition showed that they "buzzed" in a manner gov- 
erned internally by quasi-one-dimensional (plane-wave) theory, and not 
in accord with the laws of the Helmholtz ( spherical- wave) theory. A 
theory applicable to the buzzing problem was obtained by modification 
of the one-dimensional, unsteady- flow theory by introducing specified 
assumptions which greatly facilitated the computational efforts without 
impairing the usefulness of the theory for this flow. Because of the 
difficulty in computing the three-dimensional unsteady viscous flow when 
the shock was in motion ahead of the cowling and in computing when the 
shock was in motion far from equilibrium with the waves moving upstream 
from the nozzle, however, it was necessary to use a few experimental 
points at the start of the cycle as a basis on which to compute the 
many remaining pressure points of the cycle. Theoretical computations 
of the pressure-time curves, which were made on this basis for one 
supersonic-diffuser configuration with two different subsonic-duct 
lengths, agreed closely with the experimental records. The entrance 
and exit mass flow - time curves, as well as the wave diagrams, were 
also computed for these configurations by the unsteady- flow theory. 

A possible simulation of the buzzing in ram jets with combustion 
by the test of a cold-flow model in which a contraction in area replaces 
the combustion chamber is also discussed. 
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INTR0DUCTI0N 


The buzz or pulsation of axially symmetric supersonic inlets having 
central bodies has been the subject of investigations since Oswatitsch 
first observed the phenomenon in 1944 (reference 1). This buzzing is 
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essentially an oscillation of the normal shock along the spike of cer- 
tain diffusers operating at mass flows below design, and this buzz 
decreases the performance of the propulsion units appreciably (refer- 
ence 2). When a fixed-geometry ram jet is accelerated to the design 
Mach number it is necessary to fly through a range of mass flows less 
than the design value. Consequently, the importance of either avoiding 
the conditions causing the buzz or reducing and controlling the ampli- 
tude of the oscillation is evident. To accomplish this end it is, of 
course, desirable to understand the complete nature of the phenomena. 
Antonio Ferri (reference 3) has described one of the mechanisms which 
initiate the buzz. This paper continues the analysis of the buzzing 
cycle by means of a modified one— dimensional approach. The buzzing 
referred to herein is not to be confused with the pulsations of inlets 
sometimes observed due to improper operation of the combustion system. 

The causes of this buzz, in the absence of combustion, have been 
studied previously and include effects due to the presence of separation 
inside the entrance to the cowling. The choking resulting from this 
separation reduces the internal mass flow and forces the normal shock 
to move forward, so that the buzz is initiated. Reference 3 has shown 
that this separation is often due to a vortex sheet near the inner wall 
of the cowling lip. This vortex sheet arises from the intersection of 
either the conical and normal shocks or, in some cases of separation on 
the central body giving rise to lambda shocks, from the intersection of 
the legs of the lambda shock. In certain instances, the separated flow 
from the central body has appeared sufficient in itself to choke the 
cowling and precipitate the pulsations. As the buzzing cycle progresses, 
the normal shock moves rearward and assumes a position where it may dwell 
for an appreciable time, depending on internal-flow considerations, 
before the start of the next oscillation. 

Some very general theories based on internal mass- flow consider- 
ations have been expounded which give little insight into the aerodynamic 
nature of the phenomena making up the complete buzz cycle (references 3 
and 4). An attempt has also been made (reference 5) to predict the fre- 
quency only of the buzz by considering a ram jet to be similar to a 
Helmholtz resonator. The application of a quasi- one-dimensional theory 
to the buzz problem would seem proper from considerations of the ram- jet 
geometry. First of all, the area variations are gradual because of the 
necessity for high diffuser performance. Furthermore, the ratio of 
maximum area of the unit to the minimum area usually varies from 3 to 5; 
whereas for the passages of the internal-flow system the ratio of the 
length to the maximum diameter is generally in the range of 6 to 15 . 

These proportions suggest the use of an unsteady- flow theory, such as a 
quasi- one-dimensional one, which assumes a slight variation of the moving 
pressure waves from the plane form, rather than a Helmholtz type theory 
based on spherical waves (that is, large, abrupt area changes). 
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Consequently, an appropriate modified one-dimensional nonstationary- 
flow wave theory has been derived at the Langley Laboratory. This theory 
assumes quasi-steady flow for the reflections of upstream-moving waves 
which strike the normal shock and for waves passing through an area 
change. 'Because of the difficulty in computing the three-dimensional, 
nonsteady phenomenon which occurs as the normal shock moves rapidly out- • 
ward and inward during part of the cycle, it is necessary to use a few 
experimentally determined values of the waves generated by this motion 
as a basis on which to compute the remainder of the cycle. For the par- 
ticular cycles computed in this paper, four or five experimentally deter- 
mined instantaneous pressures are used to compute one hundred or more 
additional pressure points required to complete the cycle. Once the 
initial problem has been specified by the experimental points, however, 
the remainder of the solution is easily obtained from graphs with a mini- 
mum of computation and good degree of precision. 


SYMBOLS 


a 

b 



e 



h 

i 


j 

k 

m 


local speed of sound 

state in ram jet at start of buzz 

states on opposite side of typical entropy discontinuity 
effective one-dimensional values. 

states on opposite side of entropy discontinuity entering exit 
nozzle 

state ahead of typical downstream wave 

station at cowling inlet (total cross section enclosed by 
inlet ) 

state behind a typical downstream wave 

state behind reflection from a typical downstream wave 

mass flow 

mass flow at infinity through a stream, tube of diameter equal 
to the cowling-lip diameter 
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exit nozzle (measured, not effective, value) 
local static pressure 

ft 

/ 

state behind an upstream wave generated by entropy discontinuity 
passing through nozzle 

\ 

state behind second entropy discontinuity entering exit nozzle 

time, seconds 

local fluid speed 

coordinates along ram-jet axis 

area measured normal to ram- jet axis 

characteristic parameters for steady flow 


distance from cowling inlet to exit nozzle of ram jet 
Mach number 




> characteristic parameters for unsteady flow 


u 



characteristic parameter reflected at cowling by parameter 

for case when both parameters are evaluated in fluid having 
entropy equal to "effective value" . 

gas constant 

nondimens ional entropy referred to free-stream datum 
local static temperature 

coefficient of specific heat at constant pressure 
coefficient of specific heat at constant volume 
free-stream state 
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a region between flame front and nozzle of ram jet 

(B region between area contraction and nozzle of ram jet 

7 ratio of specific heats, assumed equal to 1.400 

6 distance from flame front or area contraction to exit nozzle 

angle included between ram- jet axis and line from apex of 
spike to cowling lip 

p local static density 

0 state obtained by isentropic deceleration from static condi- 

tions to zero velocity 

1 annular station at cowling inlet 

2 station immediately behind cowling in simplified ram jet 

3 station immediately ahead of nozzle in simplified ram jet 

Except for the cases noted in symbol list, subscripts denote the 
evaluation at the state or station designated by the subscript. 

THEORY 

BASIC THEORY 


One-dimensional unsteady- flow characteristic theory is governed by 
the differential equations below, which determine the variation of the 
parameters P and Q, respectively, along the corresponding character- 
istic lines having slopes of — = u + a and — = u - a 

dt dt 

5P _ a D 1 °ge A , y-l a D | | la S | (l) 

St a 0co Dt 7 a 0oo Dt 7 a Qoo St 

s s 

8 Q = a DlogeA ' 7-1 a P R,1 a & R (2) 

St a 0co Dt 7 a Qoo D t 7 a^ St 
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P s 



u 


Cfeo 


(3) 


Q = 


2 a 

y - 1 a ° M 





where + (u ± a) is the 

St ot ox 


characteristic and 
ular fluid particle 


p( ) . a( ) 

Dt 8t 


u 


3( ) 

8x 


total variation along a 
is the variation of a partic- 


These equations are derived in appendix A by combination of the 
equations of energy, momentum, continuity, and state. Equations of 
similar type may be found in references 6 and 7 in which, however, the 
area is not treated as being time dependent. 

The parameters P and Q have much the same relation to unsteady- 
flow theory as the parameters Cj and Cjj have to the Prandtl- 

Busemann characteristic theory for steady flow (for example, see refer- 
ence 8). In the steady- flow system, Cj and Cjj are invariant along 

characteristic lines of the first and second family except when the 
characteristics pass through rotational flow fields. Similarly, in one- 
dimensional unsteady flow, the values of P and Q are constant along 
the characteristics except when influenced by entropy changes or area 
changes (the latter factor influences the parameters by accounting for 
the variation of the flow from the exact one-dimensional flow). .In 
steady flow the inclination of the streamlines is proportional to the 
differences of Cj and Cjj whereas the velocity is proportional to 

their sum. Again there is a similarity in one-dimensional unsteady flow 
where the particle speed is proportional to the difference of P and Q 
whereas the local speed of sound is proportional to their sum. 


In the use of the one-dimensional characteristic equations above, a 
very lengthy and tedious step-by-step process is required to obtain a 
solution. Consequently, certain simplifications and assumptions are 
made to obtain a modified theory which appreciably facilitates the pro- 
cedure necessary to obtain a solution. 
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MODIFIED THEORY 
Basic Assumptions 


Substitution of discontinuous area changes for area gradients . - 
In order to simplify the solution, gradual variations in cross section 
are accumulated and treated as discontinuous area changes. Thus, the 
actual ram jet as shown in figure l(a) is assumed to have the shape of 
figure 1(b) where the only area changes appear at the cowling and the 
nozzle exit. The penalty for a considerable reduction in computing is 
a time distortion in the wave pattern since both the duration and time 
of generation of reflections is greatly reduced. In other words, the 
unmodified solution yields gradual reflections as soon as an area 
gradient is encountered by a wave. No reflections, however, can arise 
in the simplified ram jet until a wave hits either end of the duct, 
whereupon stronger reflections, representing the accumulation of the 
aforementioned gradual reflections, are generated. 

Restriction requiring equality of entropy of a fluid particle at 
both entrance and exit of hypothetical ram jet .- In order to simplify 
further the procedure, all waves moving through the duct proper are 
assumed to be isentropic. Thus, the entropy increase in the fluid is 
assumed to be the cumulative result of the conical shock, normal shock, 
and separation on the central body preceding the main duct. The entropy 
of entering fluid particles varies only as these factors vary and then 
this entropy (which may have a different value for various particles) 
remains constant for a particular fluid particle flowing along the duct. 

Use of quasi-steady flow to obtain reflections at duct exit arid 
entrance .- The assumption of quasi-steady flow is introduced to provide 
a simple solution for the integrated form of the characteristic equa- 
tions when waves pass through area changes or encounter the normal 
shock. This modification is based on the fact that steady flow is 
established at an area change in a very short time interval after the 
passage of the initial wave. Consequently, all the minor transmitted 
and reflected waves may be accumulated with the major transmitted and 
reflected waves, and then the group treated as consisting of only one 
transmitted and one reflected wave. Of course, this treatment also 
yields a wave pattern with a slight time distortion, but, in view of 
the labor saved, this distortion is justifiable. 

The reflection condition at the nozzle exit is then easily deter- 
mined. Since the subsequent analysis is to be limited to supersonic 
flight for which the pressure in the duct is always greater than the 
stagnation pressure required to produce sonic velocity by isentropic 
expansion to atmospheric pressure, the boundary condition of sonic 
velocity in the nozzle is dictated. Consequently, since the ratio of 
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the area of the duct to the nozzle area is invariant, there is auto- 
matically imposed the condition of constant Mach number (which must be 
equal to the effective Mach number at the start of buzz) in the ram jet 
just ahead of the nozzle under the assumption of quasi- steady flow. 
Thus, from the definitions of P and Q, one finds for the values of 
these parameters existing just ahead of the nozzle 


P 3 “ Q 3 


M^ = — = 


p 3 + q 3 


or 


( 5 ) 



5 - M e 

5 + M e 


for 7 


1.400 


(6) 


where M e is the effective Mach number of the duct. Thus, the value 
of the parameter P^ for a characteristic moving downstream at the 
nozzle uniquely specifies the parameter the reflected-upstream 

characteristic. 

For the case of reflections from the cowling inlet of upstream- 
traveling waves, the parameter P can be evaluated if the mass flow 
entering the cowling is known. During a large part of some cycles the 
normal shock is inside the cowling and therefore the mass flow may be 
determined from conical-flow theory and the external geometry. For 
that part of certain other cycles when the normal shock is ahead of the 
cowling but has not progressed out as far as the point where buzz starts 
experimental steady- flow values of the mass flow can be used. Once the 
value of the mass flow is assumed, then there exists a relation at the 
cowling entrance between P, Q, the entropy S/R, and the effective 
area A e for quasi-steady flow. As an illustration, assume that the 
external geometry is such that the conical shock intersects the cowling 
lip (that is, diffuser operating at design M^) while the normal shock 
is swallowed. Then the following equations of continuity, energy, and 
so forth, apply for quasi-steady flow: 

A-iPcAc 3 -*, = A eP 2 M 2 a 2 (7) 
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= a. 


(8a) 


7 - 


ifclf + /f 


7 - 


^ • feF • «“ 


Introduction of the entropy term into equation (7) leads to: 


7+ 1 q q 

M 2 (& 2 \7- 1 A, 


= ~L e R 

^-e 


In terms of P and Q, equations (8b) and ( 9 ) become 


(9) 


7 - 


l(^if + (-)V 2 ♦ « 2 )= 


( 10 ) 


2 _ 

|7-1 


-(t) 


7+1 
*7-1 A, 


So-S 
£1 00 


M„ie B 


(ID 


Equations (10) and (ll), plus the restriction that u > 0, which 
requires P > Q, uniquely determine the value of P for a given Q 
and S/R, or vice versa. 

Before the reflected P^' can be determined, however, the value 
of Q 2 m ust be found from the corresponding value of along the 

same upstream characteristic (fig. 2). This value of may be 

obtained from equation (2) which states, for A = constant and for 

isentropic flow of each elemental fluid volume I = 0), 


Q 3 • Q 2 = f 3 7 d | = 7 51 x 

j Jo y a o_ K y £ — a o_ K 


( 12 ) 
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Now, for the low Mach numbers encountered in a ram- jet diffuser a ^ ± 

a 0 oo ’ 

equation ( 12 ) becomes 


Similarly, 


Q. - Q 2 A 1 ^ 
3 2 7 a o„ 



AS z. 1 S 3 ~ S 2 

R 7 R 


(13) 



(14) 


Therefore, since — is known from the starting conditions, and for a 
So 

given value of P 2 and Q 2 are specified, one can determine from 

equations ( 10 ), (ll), ( 13 ), and (l4), the value of the reflected charac- 
teristic parameter Eg'. The detailed procedure to find Pg‘ for a 

specific Qg . would be the following (see fig. 2 ): 


(1) Equation (10) is solved for P 2 in terms of Q 2 . 

(2) This value of- P 2 is substituted in equation (ll) so that it 

becomes an equation of the form: Q 2 = Q^^. 


(3) This equation is then plotted together with equation ( 13 ), 
which is also in the form Q 2 = Q 2 (^), and the intersection of the two 
curves determines the correct value of S 2 /R and Q 2 . 

(4) The value of P 2 is then found from Q 2 and equation (10). 

(5) P 3 1 is finally evaluated by equation (l4). 

The most convenient method of obtaining the values of Pg» 
from Qg, however, has been found to use a plot of Pg» aginst Qg 

where Pg ' and Qg are computed by adjusting each point of the P 2 

against Q 2 fe) curve of equations ( 10 ) and (ll) by a factor of i 

\ n / 7 R 
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This manner of determining the reflection at the forward section 
of the ram jet is inapplicable when the motion of the normal shock is 
far from a quasi-steady, or near equilibrium, state with the waves which 
have come upstream from the nozzle. In other words, when the shock has 
been set in motion from a point on the external spike owing to three- 
dimensional waves emanating from the cowling and then the shock enters 
the cowling, it may not be in equilibrium with the one-dimensional waves 
existing inside. Another case where the method is inapplicable would be 
the time when separation generated at the foot of the shock causes the 
shock to move, regardless of the upstream waves from the throat. 

In addition, the use of some three-dimensional unsteady- flow char- 
acteristic theory in the presence of viscosity would be necessary to 
determine the reflection pattern from the normal shock when it is in 
motion ahead of the cowling. For these reasons, experimental pressures 
are used to prescribe the downstream waves generated when the normal 
shock is in motion either ahead of the cowling or else too rapidly for 
quasi-steady flow. In the cycles computed in this paper, both condi- 
tions occur simultaneously and necessitate the specification of only a 
few waves. 

Reference of characteristic parameters to the basic entropy and 
subsequent neglect in wave diagram of entropy discontinuities in duct .- 
Although the variation with time of the entropy of the fluid entering 
the duct is allowed, the sweeping modification is employed of referring 
all the characteristic parameters to the basic effective entropy (that 
is, effective entropy of the fluid present in the duct at start of 
buzz). This assumption drastically reduces the computational effort, 
since it eliminates the determination of an entropy-discontinuity 
streamline every time a wave strikes the normal shock and, then in turn, 
the determination of reflected and transmitted waves each time a wave 
impinges on these entropy discontinuities. It is also consistent with 
the assumption that % 1 in equations ( 13 ) and (l4) relating the 

Qoo 

increment of P or Q to the entropy variation across a discontinuity, 
since that assumption also implies the neglect of waves reflected from 
such discontinuities. 

The use of a basic entropy simply amounts to the insertion of an 
imaginary lamina of fluid having the "basic entropy" in a fluid of 
another value of entropy in order to determine the pressure and velocity, 
but not the density, of the flow. This procedure is permissible because 
the last terms in equations (l) and (2) represent the increment in the 
characteristic parameters P and Q across an entropy discontinuity 
when the pressure and velocity are constant. This fact may be proven 
as follows with reference to figure 3. By applying the definitions 
of P and Q to the characteristics in. figure 3, the following equa- 
tions are determined: 
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P d " P c “ 


OJ 


( 2 a c + 

u c ^ 

2 f 

y - 1 “<=„ 

a °oo 

V - 1 a °co 

a oJ 

’ 7-1 


(15) 


2 a d u d 


2 a c u c 


Qd " Q ° 7 - 1 a Qoo a 0oo ~ 1 a n^ a r>_. I 7-1 a, 


Ooo ^00 > 


2 ( a d ~ a c) 

°00 


( 16 ) 


since u c = u^. 


Therefore, 


AP = = — § — ( a d ~ a c) 

7-1 a_ 


2 a c / a d 


7 “ 1 a o«\ a c 


- 1 


(17) 


Introduction of entropy into equation (17) yields 

7-1 


AP = Afi = 


2 a c 


7 " 1 a o f 


7-1 AS 


( 18 ) 


2 

. fc_ 

7 - 1 AS 

7 -1 

a Ooo 

2y R 

1 AS 

for 

% 1 

7 R 


a o 


7 - 1 AS . l/7 - 1 AS\ 2 . ”| 

27 r + 2\ 27 r y • * J 


AS 

R 


since P d = P Q (19) 


( 20 ) 


Now, in figure 3(b) the true values of Q and P in the fluid of 

entropy were determined as specified by equations ( 13 ) and (l4) 

S- 

from the corresponding values in the fluid of entropy —1 by adding 

R 

S — S o 

the correction ^ d R . Consequently, when they are reduced to the 


basic entropy — ^ of the imaginary lamina inserted in the fluid, the 
R 
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1 So - S d / 

factor j (that is, 

and the original values of 


— — from equation 
7 R 

and are recovered. 



must be added. 


The relationship between the pressure and sonic speed within a 
fluid of constant entropy is: 



( 21 ) 


Since the pressure and velocity, but not density, are invariant 
across an entropy discontinuity in the flow, it is permissible to ignore 
the variation in entropy and use the values of P and Q based on the 

" S 3 

effective entropy at the start of buzz for the computation of pres- 

ft 

27 

\ ^ “ 1 

sure as = ( ■ - P ■ - + - Q A and velocity as P ~ Q in any part of 

Pe VP e + <W a 0oo 2 

the fluid regardless of entropy. The wave diagram will then be greatly 
simplified, but will not be a true picture of the flow in regard to 
stream density or mass flow. 

Of course, if the entropy at a point is- known, the mass flow can 
be determined from the values of P and Q in the imaginary lamina of 
basic entropy. When reference is made to figure 3(b) where condi- 
tions c correspond to the imaginary lamina. 


^c = Pc^c 
m d p d u d 

1 7-1 s d"S c 

- (I cV . 7 R u c 

\ P d) , u d 


( 22 ) 


(23) 


S d- S c 7-1 

7 ^nce p c = p d , u c = u d 


(24) 
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When equation (24) is expanded and higher order terms are 
neglected, 


“d _ 


m„ 


= e 


s c" s d 7-1 

R 7 


* i + 21^1 !s_- S d + 
7 R 


( 25 ) 


The ratio of 



(26) 


since 



( 27 ) 


Then equation (26) becomes 


me = Ac + Q c\ y V pc - 

*b \ P b + %) \ p b - %) 


( 28 ) 


Thus the ratio of the actual mass flow to the mass flow at the start of 
buzz becomes: 


2 

\ _ “d m c Ac + QcY" 1 /^ - <*c\/ , 7 - 1 S e - S d\ 

m b m c m b'K + Qj \p b “V\ y / 


( 29 ) 


There still remains the problem of determining the local entropy. This 
task, if it is to be precisely accomplished, requires the plotting of 
the entropy-discontinuity space-time histories. Without these histories, 
however, it is possible to obtain easily a so-called "extrapolated value" 
of the mass flow entering the cowling by assuming quasi- steady flow from 
infinity to the cowling lip. This value is not a true value when the 


CONFIDENTIAL 



NACA RM L52A18 


CONFIDENTIAL 


15 


shock is inside the cowling and out of equilibrium with the upstream 
Q waves, since the mass flow entering must be equal to when the 

conical shock is at the lip of the cowling. This error is due to the 
fact that the parameters P and Q were determined experimentally at 
stations to the rear of the shock and included the effect of the shock 
motion. Consequently, use of these P and Q terms evaluated 

a t ^ = 0 give the mass flow just behind the shock instead of at the 

cowling entrance so that the values are only correct when the 'shock is 
ahead of the cowling or when it is in equilibrium inside the cowling. 
Since these parameters determine the local pressure and velocity, the 
assumption of quasi-steady flow permits the determination of the local 
temperature from the energy equation. Then, when the pressure and the 
temperature are known, the local entropy and, in turn, the mass flow 
may be found. 

When the presence of entropy discontinuities in the flow is dis- 
regarded in the wave diagram, the waves originating from the arrival of 
these discontinuities at changes in cross section are, of course, 
omitted. In particular, the upstream wave generated by the passage of 
a discontinuity through the choked exit nozzle is neglected and this 
omission must be justified. Figure 4 shows the wave pattern in the 
diffuser created by the passage of an entropy change through the nozzle. 
The following relations then apply: 



(30) 


Q g = Qp + — where — 
g i 7 r r 



( 31 ) 


p - p . as 

p g f 7R 


(32) 
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2 


7-1 


2 


7-1 


- M c 


+ M c 





= Qf + 



AS 
7 R 


( 33 ) 


( 34 ) 




Qq - Qg 


2M e S g ~ S f 
— 2 — + M 7R 

7-1 


( 35 ) 


( 36 ) 


Thus, waves of strength Qq - Qg are omitted in a solution neg- 
lecting entropy discontinuities ^ which is equivalent to the assumption 

. In cases where this assumption is not justified. 



hut another entropy discontinuity 



having opposite sign so that 


^r 

R 


- t « i 


soon follows the initial discontinuity through the area 


change, then the assumption has only introduced a time distortion in the 
wave pattern by treating a positive wave followed by an approximately 
equal negative wave as a single wave of zero strength. 


Effective values employed in the modified theory .- One-dimensional 
flow presumes no variation of the flow properties in a direction normal 
to the stream velocity. Consequently, in order to apply such a theory 
to an actual ram jet where boundary- layer separation and centrifugal 
force actually cause a variation in the fluid properties across a sec- 
tion of the flow, effective values of the fluid properties and cross- 
sectional area must be employed. The effective values are defined as 
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the values satisfying the one-dimensional channel-flow relationships 
regarding continuity, energy, momentum, and the equation of state. 

Since the effective area and Mach number A e and Mg have such 

an important role in the theory, the equations which determine them are 
presented below. Steady flow is assumed to exist for an infinitesimal 
time before the buzz is started. 


Mass flow at start of buzz = 11 ^ = p e M e agA e = P^M^a^A^ 


^b 

®00 


(37) 


M =M 

” e ^ a e P e A e m* 


(38) 


From the definition of entropy and equation of state. 


2 S e" S <» 


- ST 1 


(39) 


Since a n = a n at start, equation (39) becomes 

U OO 



Combining equations ( 38 ) and (40), 

7+1 
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Insertion of energy equation into equation 


(4l) gives 



M 


e 





(43) 


The specification of Ae plus experimental values for the other terms 
on the right-hand side of equation (43) prescribes M e uniquely. 

Approximation for simplifying construction of characteristic 
diagram . - Instead of the step-by-step procedure of constructing the 
characteristic lines with local slopes dependent on the sum or differ- 
ence of the local fluid speed and sonic speed, the value of these speeds 
at the ends of the duct is used to obtain an average slope. Thus char- 
acteristic lines of the first family are drawn with a slope 



(44) 


evaluated as the average of the velocities for — = 0 

L 

Similarly, for the second family, 


and 


x 

L 


= 1 . 



( 45 ) 
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whereas for particle space-time histories 



( 46 ) 


A linear variation of the parameters is used for interpolation 
between characteristics of the same family. 

The characteristic theory has now been so greatly simplified that 
the solution is easily obtained with the help of three plots, namely, 

2 

(a) pressure-sonic velocity, p = p^JL ^ 7 " 1 (b) P - Q reflections at 

cowling inlet and (c) Q - P reflections at the sonic nozzle. After 
the experimental values of the few waves generated by the external 
motion of the normal shock are plotted, the remainder of the solution 
consists of reading values off the curves and plotting the 
characteristics. 


APPARATUS 


The apparatus employed for the experiments is shown in figure 5 
and is essentially the same as that used in reference 3 with certain 
modifications necessary for the conversion from steady- flow to unsteady- 
flow measurements. The total-pressure rake at the survey plane 
(x = 11.50 inches), figure 5 , was removed and an instantaneous static- 
pressure pickup employed in its place. The instantaneous pressure 
pickups were MCA miniature electrical pressure gages mounted flush 
with the diffuser wall. Also, since only reflections from the nozzle 
exit were desired, the mass-flow orifice was eliminated. 

Tests were made at a Mach number of 1.9 in the Langley 4- inch by 
5 -inch blowdown open jet by using low-humidity air from pressurized 
tanks. The test Reynolds number based on cowling- inlet diameter was in 

ft 

the range between 2.9 x 10 u and 3.6 x 10 . Three different subsonic- 
duct lengths of 34 inches, 156 inches, and 390 inches were employed with 
instantaneous pressure pickups located. along the diffuser as shown 
schematically in figure 6 . The 90° elbows necessary to attain the 
desired duct lengths within the space limitations of the laboratory are 
indicated. The length of the central bodies (which had the same after- 
section in all cases) relative to diffuser lengths is also shown. 
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Tests were made on the three supersonic-diffuser configurations 
listed in the table below: 


Config- 

uration 

Semicone 

angle 

(deg) 

Maximum diameter 
of central body 
(in. ) 

C owling 
(reference 3) 

Values of cowling- 
position parameters 

N 

I 

25 

0.80 

E 

44° O’, 40° 40', 36° 50’ 

II 

25 

1.10 

D 

43° 10*, 39° 50*, 36° 50* 

III 

30 

0.80 

E 

48° 30*, 43° 40 *, 40° 4o* 


The dimensions of these central bodies and cowlings are given in 
reference 3 with the exception of the central body of configuration I, 
the dimensions of which are given in figure 7* 

A commercial concentrated arc light served as the source for a 
shadowgraph system in which the image on a ground-glass screen was 
recorded on either of two 16-millimeter motion-picture cameras. For 
the longer diffuser lengths a low-speed camera operating at 64 frames 
per second was used, whereas for the short diffuser a Wo’llensak Optical 
Company Fastax camera running approximately 500 frames per second was 
employed. Signals from the pressure gages were amplified by a 
Consolidated type 1-113 amplifier and then recorded on sensitized paper 
in a Consolidated type 5-114 recorder employing galvanometer elements. 
The frequency response of this system was flat from 2 to 300 cycles per 
second. In order to correlate movie film with the pressure records, a 
timing light was attached to one edge of the shadowgraph glass screen 
and the circuit energizing this light connected in parallel with one of 
the recording galvanometers so that every time the light flashed a 
"blip" appeared on the pressure record. A synchronizing device was 
attached to the low-speed camera which flashed the light on the shadow- 
graph screen once each frame, while 60-cycle alternating current was 
used to operate the light when the Fastax camera was used. In order to 
reduce the labor of correlation, an interrupter was placed in the light 
circuit to stop the flow of current for a noticeable period four or five 
times a second, and hence provide blank spaces on the film and pressure 
record which could be easily counted. 

A manually operated, 2-inch gate valve served as the exit throttle. 
Throttle position was recorded on the pressure film by means of an elec- 
trical circuit employing a ten-turn linear' potentiometer attached to an 
extension of the valve stem. 
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TEST PROCEDURE 


Tests of the various configurations were made by the following pro- 
cedure: the desired central-body projection was set by scale measure- 

ment and, with no tunnel air flowing, a trace was then recorded on the 
pressure film; in the cases of runs made with the 64- frames-per- second 
movie camera, a strip of "no flow" motion-picture film was also made. 

The exit throttle was opened to such an extent that the mass flow would 
be limited only by the supersonic-diffuser configuration. Then the 
tunnel was started and brought up to operating conditions. The movie 
camera and pressure recorder were started and the throttle closed very 
slowly until buzz began, at which point the throttle position was held 
approximately constant for a few cycles before continuing on to a lower 
valve opening. This procedure of a slow throttling interrupted by 
pauses at constant valve area was continued for the duration of the 
test run except for a rapid opening of the throttle at the end of the 
run. The tank stagnation pressure as well as the differential-gage 
reference pressure were recorded during the run. A no- flow pressure 
trace was also made after a run when possible. The differential gages 
and valve-position recorder were calibrated statically each day on which 
tests were run. 

The cowling-position parameter, or angle between the axis and a 
ray from the cowling lip to the apex of the central body, was determined 
from enlargements of the no- flow motion pictures. The outlines of the 
conical body were extended to intersect at the apex, and then angles 
were measured from this point to the outline of the cowling lip. This 
method, however, has a low accuracy (of the order of 20 to 30 minutes) 
because of the presence of a slight fuzziness on the shadowgraph 
outlines . 

The value of A n , area of the valve opening normal to the axis, 

was found for each position of the throttle wheel by integrating the 
shadowgraphs obtained by placing the valve in the optic system with the 
flow axis alined with the optic axis. The value of A n is only a rela- 
tive, and by no means exact, measure of the effective throttle area 
since the valve seat is appreciably recessed. Consequently, the effec- 
tive area is always much larger than A n , but the variations in A n are 

still indicative of variations in the effective throttle area. 
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RESULTS AND DISCUSSION 
COMPARISON OF THEORETICAL AND EXPERIMENTAL RESULTS 


Theoretical computations of buzzing cycles for configuration I 
= 44° + 20‘) with ‘duct lengths of 156 inches and 34 inches have been 

made. These cycles were typical of the configurations which buzzed with 
the vortex at the lip, and, in addition, were the most complicated in 
the sense that the one-dimens iondl wave pattern included a large number 
of generated and reflected waves. The cycle for the 156-inch duct is 
discussed first because it more clearly shows the waves a nd their 
reflections since the duct is long enough to permit the complete gener- 
ation of a wave before its own reflection from the nozzle exit returns 
to alter the generation. 

The basic computations for this cycle are given in appendix B. The 
characteristic or wave diagram, including one space-time history for a 
fluid particle which enters as the buzz starts, is shown in figure 8, 
whereas the theoretical pressures are plotted on the experimental pres- 
sure record in figure 9* Flagged symbols are employed in this paper to 
denote those values computed from the experimental records. A compari- 
son of the many additional theoretically computed points with the experi- 
mental pressure record shows generally good agreement for all five pres- 
sure traces and thus offers strong evidence . that the modified theory is 
valid. 

The lowest trace of figure 9 is that of the valve-position indi- 
cator. This is not a smooth curve because of the electrical "pickup" 
and, consequently, the value employed to determine A n is the average 

of the maximum and minimum of that part of the trace where the "pickup” 
is smallest. The topmost blips indicate the timing-light flashes used 
to synchronize the motion pictures to the pressure record. 

A strip of movie film (fig. 10) the frames of which are numbered to 
correspond to the numbers in the blips at the top of the pressure record 
(fig* 9) shows the position of the normal shock at various times. The 
shock moves outward between the first and second frames and then returns 
into the cowling between the third and fourth frames, where it remains 
until the fourteenth frame. Three-dimensional optical effects are 
responsible for the apparent presence of a normal shock ahead of the 
cowling in frames 4 to 13. It is evident that the flow pattern is not 
precisely symmetrical, a fact caused by a slight eccentricity in the 
location of the central body. 


CONFIDENTIAL 



NACA EM L52A18 


CONFIDENTIAL 


23 


At a given value of x/L (in fig. 8) the increase of the charac- 
teristic parameters P or Q with time indicates the passage of a 
downstream or upstream compression wave, respectively. S imi larly, a 
decrease in P or Q denotes an expansion wave. 

The wave pattern of the cycle is now discernible from figures 8 
to 10. The normal shock moves out because of the separation in duced by 
the vortex and assumes a steady location on the spike which is dictated 
by the choking in the inlet. This motion of the shock causes an expan- 
sion wave to move downstream to reduce the mass flow. At the nozzle 
exit the downstream expansion is reflected as another expansion which 
moves upstream until it overtakes the normal shock whereupon, after a 
slight delay, this shock retreats into the cowling. The rearward motion 
of the normal shock creates a strong downstream compression wave which 
is closely followed by a much weaker expansion wave. This weaker expan- 
sion wave is probably due to the fact that when the shock moved into the 
inlet it went too far downstream. Consequently, when separation occurred 
behind it, the shock had to move up the channel somewhat and in so doing 
generated a weak expansion. In general terms, the shock overshot the 
equilibrium position when it moved back and, hence, had to adjust its 
position by moving upstream a little. It is important to note that the 
vortex from the intersection of the conical and normal shocks had to 
touch the inner wall of the cowling as the shock returned into the dif- 
fuser and yet separation did not occur to initiate the next cycle. An 
explanation for this phenomenon may be found in the following consider- 
ations. First, the normal shock is retreating so that the total pressure 
losses across it are reduced. It is then possible that inside the 
cowling the static pressure of the fluid on the outer side of the vortex 
sheet will not be required to equal its stagnation pressure; therefore, 
separation may be avoided (see reference 3). Second, the strong upstream 
expansion tends to prohibit separation by accelerating the flow and pro- 
ducing a favorable pressure gradient along the diffuser. 

The combination of a strong compression followed by a weaker expan- 
sion is reflected as a unit, with the same sign both at the nozzle exit 
and at the shock inside the cowling. Each reflection, however, lessens 
the magnitude of the waves so that, as the pressure in the diffuser is 
increased by the waves, the steps due to the waves approach zero. The 
shape of the pressure-time curve so generated has often been described 
as a "relaxation type" but is, in reality, a series of step functions 
separated by approximately constant pressure plateaus. The compression 
waves eventually force the normal shock to the cowling lip where the 
vortex attains a position from which it can start the cycle again. In 
this particular cycle, the normal shock starts to move outside the 
cowling in the beginning of the next oscillation just prior to the 
arrival of the seventh upstream reflection of the compression wave 
which by then has degenerated from an original calculated magnitude of 
4.7 pounds per square inch to only 0.1 pound per square inch. A 
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possible explanation for the origin of the second cycle between the two 
compression reflections might be that the earlier of the two waves placed 
the vortex from the normal shock in a position where a random motion 
could enable the vortex to cause separation. Another possibility might 
be a time lag required for the separation to grow to a value where it 
could choke the flow. . 


The magnitude of the errors which can be introduced by the various 
assumptions may also be determined. By using the exact formula for the 
variation of the characteristic parameters across a constant-pressure 
entropy discontinuity with the lowest value of a/a 0oo occurring in the 

cycle, the error in the reflected value of P may be found from equa- 
tions (10), (ll), and the following equations (47) and (48) (which are 
similar in form to equation (l8)): 


Q 2 - Q e 


/ Sg-Sg y_i 

2 a e 1 R 27 

7 - 1 a D ' 

U 00 



(47) 



s e" s 2 7-l\ 
. R 27 1 


(48) 


This error in P is equal to 0.0019, which corresponds to a posi- 
tive pressure deviation of approximately 0.1 pound per square inch. 

AS 7-1 

Consequently, the assumption of -- fa 1 and e ^ ^ 1 + ^ ^ 

a °«> R 27 

introduces an error no larger than 0.1 pound per square inch at each 
reflection and usually is considerably less since the low temperature 
is present for only a. small part of the cycle. 

The highest value of the entropy in the cycle is 0.2723; therefore, 
the maximum error involved by neglecting the passage of entropy discon- 
tinuities through the throat must not' be greater than the value computed 
from equation (36), that is, - Q g = -O.OO59. 

This neglected increment is approximately equal to an expansion 
wave of 0.33 pound per square inch. The entering entropy, however, 
drops enough in 0.015 second after the maximum so that the net. value 
of the reflected waves then would be less than 0.1 pound per square 
inch. Thus, negligence of entropy has, introduced a time distortion in . 
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the wave pattern by omitting an expansion wave of 0.3 pound per square 
inch followed 0.015 second later by a compression wave of 0.2 pound per 
square inch, both of which will damp out as they reflect back and forth 
in the diffuser. 

Since the strongest single compression in the cycle has a pressure 
77. ^ 

ratio of — = 1.065 (see gage 4 of fig. S a.t t w 0.04), the assump- 

72.6 

tion of isentropic waves is valid beyond question. 

The assumption of quasi-steady flow ahead of the cowling for pur- 
poses of entropy and mass-flow computations neglects the velocity of the 
normal shock. From the movie film it' appears that the shock velocity 
is about 4 feet per second, which, computed on a normal-shock basis for 
a free-stream Mach number of 1.90, yields a ratio of stagnation velocity 
of sound behind the shock to that of the free-stream stagnation sonic 
velocity of 1.00046 as the shock moves outward. Under these conditions, 
the true local entropy at the lip is greater than the steady- flow value 

AS 

by an increment of = 0.003, and the correct mass flow is obtained 

by multiplying the steady- flow value by a factor of 0.999* The possi- 
bility of other waves between the cowling and shock must also be con- 
sidered, but, since the time interval for flow to pass from the spike 
tip to the cowling is approximately 0.0001 second and none of the 
unsteady waves has a very large pressure slope dp/dt, these waves will 
also have negligible effect. 

The influence of a variation in A e 
changes M e , is shown by figure 11. The 

inches used for the sample computations ( and comparison of theory and 
experiment was the channel area normal to the axis at x = 12.75 inches. 
This area is not the largest area of the diffuser, since at the end of 
the central-body support strut the area increases to 3*339 square 
inches. Thus, an effective area of 3*190 square inches presumes an 
average boundary- layer correction of 0.1^9 square inch in the 
remainder of the diffuser. • Since there was no 'accurate method of meas- 
uring this correction, the- influence of arbitrarily varying the effec- 

• ■ * ■ ■ ' in y. 

tive area to 2.99 square inches was investigated. For a — 2. = 1.000 

m oo 

and A e = 2.99 square inches, M e was found to be O.I63. In a similar 
manner as set forth in appendix B for . A e = 3.190 square inches, the 

wave diagram -and pressure-time curves/ were, .'determined by using the 
identical few experimental pressures.-. .An, examination of figure 11 shows 
that the only significant difference vis: tha : t, the pressure .in the middle 


or m^m^, which, of course, 
A e value of 3*190 square 
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square inches than for 

1 . 000 . 


flow at start of buzz from 

m h m Vi 

— = 1.000 to O .985 for an effective area A = 3 .I 9 O square 

inches was determined. This computation was complicated by the fact 
that it was necessary to shift from the reflection curve at the cowling 

mp mo nk, 

lip for -^ = 1.000 to -&=-£= O .985 in order to complete the 

AU 00 1^00 m oO 

cycle. These reflection curves are plotted in figure 12 , along with a 
line between them which was arbitrarily assumed to represent conditions 
as the shock moves from the cowling lip out to where the buzz starts 
again at a mass flow of O. 985 . In figure 11 the untailed diamonds 

represent the pressures if the ratio ^ = 1.000 is assumed once the 

m 00 

shock has returned into the cowling until the next buzz starts. The 
tailed diamonds represent the adjustment obtained by following the 
path 1, 2, 3 (fig. 12) along the reflection curve at the cowling as the 
shock moves out. Higher pressures are evident throughout after the 
experimentally adjusted pressure portion of the cycle, with the cycle 

mp 

f° r = - 1 -* 000 resulting in greater pressures than those in the cycle 
which allows for the mass-flow reduction from ^ = 1.000 to ^ = 0.985 

m oo mo, 7 ^ 

A close inspection of figures 9 and 11 reveals the fact that the 

mu 

theoretical pressures for - 2 - = 1.000 do not return to the original 

m oo 

trip value of 87. 0 pounds per square inch, but reach a value of only 
86.7 pounds per square inch. This discrepancy probably is due to the 
continuous slow throttling of the exit valve, since the valve— position 
indicator showed a 3 . 4-percent decrease in A n during the cycle from 
its value at the start of O .67 square inch. Such a throttling could 
generate upstream compression waves to increase the experimental pres- 
sure above the theoretical value computed for constant throttle area 
and Mg. 

The experimental pressure record also shows a pressure loss through 
the pipe plus an increasing hash or noise level as x/L increased. 

These results were to be expected since there were three 90° elbows in 
the standard pipe employed as a duct and in addition a chatter in the 
gate valve which was used as a throttle. The entropy increase along 
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part of the cycle is higher for Ae = 2.99 
A„ = 3.19 square inches and the same — = 

m oo 

Also, the effect of changing the mass 
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the pipe, as well as the waves introduced by the chattering valve, of 
course, were neglected in theoretical work and are, therefore, another 
source of error. 

Approximate values of the mass flow entering and leaving the ram 
jet are presented in figure 13. The "entering" mass-flow curve is 
really an extrapolated curve determined from pressures downstream which 
gives the mass flow behind the normal shock, and consequently, when the 
normal shock is in motion inside the cowling, an error will be present. 
This error is evident in the apparent decrease in entering mass flow 
from t = O.O36 to t = 0.0375 which, in reality, is just the shock 
moving upstream inside the cowling after overshooting when it returned 
from the external spike. The mass-flow curve is therefore adjusted as 
shown by the dashed line to give the correct entering mass flow. 

An approximation for the exit mass flow was obtained by computing 
the mass flow on an effective entropy basis (plain squares in fig. 13) 
and then correcting this value when the fluid which entered the cowling 
at the start of buzz had reached the exit. Each value of the mass flow 
at the exit at time t^ was multiplied by a correction factor based on 

the entropy existing at the cowling inlet at time tn - At, where At 

was the time required for fluid to flow from the cowling to the exit at 
a velocity determined from the effective Mach number Mg and steady 
flow. 


At = 


— 13^0 _ 0.0763 sec 

0.1528 x O.99767 x 1117.4 


M, 


a. 


(49) 


Integration of the area under the mass-flow curves and division by 
the product of mass flow at start of buzz multiplied by the cyclic time 
yields the average mass-flow ratio: 


f m(t) dt 


jn 

m v 


v 


(50) 


A value of 0.951 is obtained for the adjusted curve for the cowling 
inlet and 0.953 for the nozzle exit (if a time lag At* = O.O83 second 
which is the time required for the initial fluid to reach x the exit, is 
used, a value of O.952 results at the exit). According to these approxi- 
mations, which include the mass adjustment shown by the shaded area in 
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figure 13, there is slightly less mass in the diffuser at the time of 
the start of the second cycle. This condition is a logical one, since 
air of higher entropy than the effective value is still flowing in the 
diffuser at time the second buzz starts. The assumptions made to com- 
pute the mass flow, however, are so rough that the accuracy of the third 
significant figure is very doubtful. 

The same central body and cowling were tested with the short duct 

which prescribed an over-all length of 3^ inches, or ratio of length to 

maximum internal diameter of 16. 5 approximating the proportions of long 

ram jets. The wave diagram and the pressure record containing both 

theoretical and experimental points are shown in figures 14 and 15 for 

the fourth oscillation after starting. This oscillation was chosen 

because of its closer similarity to the one for the longer length in 

order to facilitate comparison. Since the value of p for this run 

o 00 

was 100.7 pounds per square inch absolute, compared to 101.0 pounds per 
square inch absolute for the longer diffuser, the value of p^ was 

reduced from 87. 0 pounds per square inch absolute to 86.6 pounds per 
square inch absolute. The theoretical computations are based on 

A e = 3-190 square inches, ^ = 0.993, M e = 0.152, but the cycle is 

m 00 

not, however, adjusted for the motion of the shock from = 1.000 

moo 

IUp 

to — £• = 0.993 just prior to the next buzz. The value M„ = 0.152 
00 

mu 

was chosen after an initial computation employing — = 1.000 and 

m oo 

M e = 0.153 yielded results consistently low by 0.5 pound per square 
inch in the latter half of the cycle. The change in the cowling-lip 

mu im 

parameter to reduce the mass flow from — = 1.000 to — = 0.993 is 

^ co moo 

less than 0° 10’ (see fig. 25, reference 3) which is much smaller than 
the measurable accuracy of the parameter by the methods employed. Con- 
sequently, such an adjustment does not invalidate the solution in any 
way. 


Again, close agreement between theory and experiment is attained. 

It is evident that the compression waves computed lead the experimental 
ones by about 0.001 second. This error is due to the initial time 
location of the first adjusted compression wave being in error by that 
amount, but, since the response of the recording apparatus was flat to 
only 300 cycles per second, this error was not worth correcting. The 
valve chatter is much more pronounced in the third trace where the pres- 
sure pickup was located only 3^- inches' ahead of the valve. 
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Points indicating the arrival of similar waves at gage 1 on the 
first trace of figures 9 and 15 are numbered to show the correspondence 
of the wave patterns for the long and short diffusers. The magnitude of 
the pressure variations is smaller for the short configuration because 
the reflection of the expansion wave, generated when the shock moves out 
to start the buzz, returns to hit the shock and stop its motion before 
it can move out far enough to reach a new equilibrium position on the 
spike. From figure 9 it is evident that the shock will move out for 
about 0.013 second before reaching equilibrium if its outward motion 
is due to vortex- induced separation and it is not forced out further by- 
being hit by a strong compression. In contrast it takes only about 
0.005 second for a wave to traverse the short duct in both directions. 
Consequently, for the same value of A Q , all waves have less strength 

for the short diffuser than in the long one. The maximum and minimum 
pressures for the cycles computed were 87. 0 pounds per square inch and 
72.6 pounds per square inch in the long configuration contrasted to 86.6 
and 79.O in the short one. Similarly, the minimum entering mass flow, 
referred to free-stream cowling area, was 51 percent and 79 percent for 
the long and short diffusers, respectively. 

The moving-picture film (fig. l6) taken concurrently with the pres- 
sure record and numbered to correspond to the numbers on the pressure 
record shows the shock motion to be much smaller for the short pipe th an 
in the long one. The shock motion also agrees with the theory in every 
respect. 

The curves of mass flow against time (fig. 17) were computed by the 
same method employed previously. The ratio of average mass flow to mass 
flow at start of buzz is 0.979 at entrance and 0.977 at exit. The 
entering value was not adjusted for the motion of the shock from 
m ? m h 

— = 1.000 to ^ = 0.993 prior to the next cycle which would reduce 

m oo m OO 

the average entering value slightly. The ratios of the average mass 

flows to the mass flow based on 100 percent capture area — are 0.972 

moo 

and 0.970, respectively, compared to 0.951 and 0.953 for the longer 
diffuser. 

Since all the variations in p, a, S/R, and so forth are smaller 
for the shorter pipe, the errors introduced by the various assumptions 
(except the one made by concentrating the area changes) which simplify 
the exact one-dimensional theory are less than for the longer pipe. 

The distortion in the wave pattern due to assuming that all the area 
changes are concentrated at the inlet and exit is increased in the 
smaller diffuser since the actual area variations occupy a larger 
percentage of the- total length (see fig. 6). 
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GENERAL NATURE OF THE BUZZ CYCLES 


The manner in which the cyclic period is usually changed due to 
throttle adjustment is shown in figure 18, which is another part of the 
pressure record for configuration 1 , 9 ^ = 44°, L = 156 inches. At 

point I (see fig. 18 (a)) the "basic cycle" is (l) a downstream expan- 
sion, (2) an expansion reflected upstream from the exit which pulls the 
normal shock back in, (3) a downstream compression, and (4), a reflected 
upstream compression which forces the shock out to start the next cycle. 
Two smaller compressions are superimposed on the basic cycle, since the 
strong upstream compression which starts the next cycle passes through 
the contraction at the cowling entrance and reflects a weaker compression 
which shows up in points 1, 2, and 3- This is, in turn, reflected and 
re-reflected as a compression upstream from the exit, points 4 and 5, 
downstream from the cowling contraction, points 5 to 7> and finally 
upstream again, points 8 to 10. Further reflections are negligible. 

As the throttle is opened, the reflections from the throttle become 
smaller. This phenomenon is due to that fact that, if the conditions h, 
j, k represent, respectively, the conditions ahead of a downstream 
wave, behind the downstream wave, and after the upstream reflection of 
the wave from a sonic valve, the ratio of the wave strengths for isen- 
tropic waves is by linearization approximately: 



(51) 


As the sonic area is increased, the Mach number must increase since 

the ratio of the sonic throat area to A^ becomes larger and the 


reflected pressure wave becomes weaker approaches 1.0 1. 


If each reflection becomes weaker while the strength of the initial 
expansion remains approximately constant, it will take an increasing 
number of compression reflections to move the normal shock to a point 
where the vortex starts the next cycle. Figure 18(b) shows this clearly 
where I requires only one compression reflection; II, two; III, three; 
IV, four; and V, five. In some cases it is evident that the cycle starts 
between reflections and, as explained earlier, starting at this time is 
probably due to the previous wave moving the normal shock to a position 
where the vortex is on the verge of causing separation so that any random 
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variation in the flow starts the cycle. Note that the basic cycle is 
approximately the same for all the frequencies and that the minor com- 
pression wave becomes unnoticeable as the valve is opened. 

The pressure record obtained when the valve was quickly opened for 
the short diffuser (fig. 19) has similar characteristics. The waves, 
however, are not as sharply defined as in the longer pipe. 

The short duct shows a marked decrease in the wave amplitude as 
the valve was opened. The magnitude of the upstream compression which 
knocks the shock out of the cowling determines the expansion strength, 
since a strong compression will knock the shock further out before it 
can be stopped by the reflected expansion. In the longer pipe, however, 
(fig. 18) the amplitude varied only slightly, since the normal shock 
moves out until it reaches an equilibrium with the sepa rat ion- choked 
cowling before any of the reflected expansion returns to alter its 
motion. , 

An over-all examination of the results of the nine different 
supersonic-inlet geometries tested revealed a general behavior pattern 
for the buzzing after it was once initiated. This result was time 
although the 30° configurations were always choked at the cowling inlet 
so that the normal shock could not enter during steady flow. Also, con- 
figuration II for = 39° 50’ buzzed due to the vortex from the 

separation- induced lambda shock striking the cowling. 

The pressure records for configuration III were very similar to 
those shown for configuration I = 44° after the cycle in which the 

shock motion first be'came appreciable. This cycle was usually preceded 
by a few short oscillations of much lower magnitude. 

The pressure records for configuration II were similar to those for 
configurations I and III only after the amplitude of the buzz had become 
large enough for the basic cycle to manifest itself. With this body 
and 0j = 43° 10’, the buzz did not start until the normal shock was 

appreciably more than halfway out on the spike and the vortex appreciably 
inside the cowling. The buzz then began in a manner best described as an 
unstable organ-pipe type of resonance. As evident from figure 20, the 
oscillations in the short diffuser increased in amplitude rapidly and 
then remained constant. Moving pictures showed that the buzz was not 
started by any vortex action, but rather by a boundary- layer instability, 
and that the magnitude of the oscillations built up until the normal 
shock was moving from a point one-quarter of the way down the spike to a 
point off the spike ahead of the tip. Shadowgraph excerpts (the frames 
are numbered to correspond to the numbers in fig. 20) of this configu- 
ration are shown in figure 21. The table on the following page gives 
the difference in pressures from the peak to the valley at each of the 
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three pressure pickups at times I, II, and III for three consecutive 
oscillations noted in figure 20: 


Time 

Ap for | 

(psi) 

gage 1, ^ = 0.338 

gage 2, £ = 0.721 

J_i 

gage 3, = 0.897 

I 

7.7 

10.7 

12.5 

II 

9*5 

13.1 

14.7 

III 

11.0 

15.6 

18.2 


This table definitely refutes the use of a Helmholtz resonator 
analysis even for this cycle, which in the absence of starting due to 
the vortex is nearest to a sinusoidal cycle, and hence would be most 
favorable to Helmholtz theory. At the same time, it substantiates the 
plane-wave theory for the short diffuser having a fineness ratio of 
less than 17. The basis for this conclusion may be found in the primary 
assumptions employed in the Helmholtz derivation; namely, the pressure 
is uniform at any time in the large cavity, and thus the pressure fluc- 
tuations must have the same magnitude throughout the cavity. The pres- 
sure records, however, show that in the duct, which in the Helmholtz 
analysis (reference 5) is considered the resonator cavity, the magnitude 
of the fluctuations varies by about 50 percent in the same cycle. Fur- 
thermore, although the oscillation is unstable and is increasing rapidly, 

the pressure variation near the valve = O.897) during one cycle is 

still larger than the pressure variation further upstream = 0. 



during the next cycle. In other words, the diffuser behaves similarly, 
to an organ pipe with relative nodes in velocity and anti-nodes in pres- 
sure at the exit for this cycle. Additional verification was obtained 
from the longer version of the same supersonic-diffuser configuration 
where a typical oscillation generated a pressure variation of 1.2 pounds 

per square inch at -•= 0.074; whereas at - = O.705 the variation was 

L L ■ 

4.5 pounds per square inch. 


Figure 20 also is an exception to the general . trend of increasing 
frequency with decreasing valve area. since the periodic time increases 
from 0.0149 second to 0.0155 second during resonance. The reason for 
this result is found in the "inertia" effects of the external motion of 
the shock. From figures 14 and.l6 it is obvious that the shock does' 
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not start to return to the cowling as soon as struck by the upstream 
expansion wave (t = 0.0114 second, frame 5)> but delays before starting 
inward (t = 0.0152 second, frame 7). As the amplitude of the shock 
motion increases, more of the expansion will have to strike it to 
initiate its movement toward the cowling. Thus, for figure 20, the 
periodic time, comprised of the lag time and basic cycle time (that is, 
that time required for an expansion and compression to traverse the duct 
in both directions), increases as the amplitude and lag increase. 

It is then obvious that, for any finite-amplitude oscillation gov- 
erned by reflections from the exit, the periodic- time must be greater 
than the basic cycle interval. The minimum period for the configura- 
tions tested was 0.0148 second, with a basic cycle time of 0.0101 sec- 
ond, when waves reflected from the exit controlled the buzz. Configu- 
ration II, 0^ = 36 ° 50’, however, oscillated at a very high frequency 

with a period of approximately 0.002 second when the exit area was 
reduced appreciably. This oscillation may have been governed by the 
supersonic diffuser entirely, and perhaps resulted from the periodic 
boundary- layer variations passing through the first minimum. 


A POSSIBLE APPLICATION TO SIMULATE CONSTANT-AREA HEAT ADDITION 
BY COLD FLOW WITH AREA CONTRACTION 


Once the quasl-one-dimensional nature of the wave phenomena has 
been established, the possibility arises of employing area changes to 
simulate heat addition. An examination of equations (l) and (2) shows 
that the parameters P and Q are functions only of the changes ini A 
or S/R. It is possible therefore to get the same variation for one 
(not both simultaneously) of these parameters due to a change in S/R 
while A is held constant that can be obtained by holding S/R constant 
and varying A. A generalization of this principle is used to produce 
theoretically a similar wave pattern for the case of constant-area heat 
addition and isentropic area contraction for configuration pairs like 
those shown in figure 22. 

Consider the simplified combustion process in a semi-infinite pipe 
diagrammed in figure 22(a) where the combustion is assumed to occur in 
a discontinuous flame front. If conditions at the flame front may be 
assumed to reach a steady or equilibrium flow after each transmission 
and reflection of a pressure wave, the. equations of continuity, energy, 
momentum, and state may be incorporated in a trial-and-error process to 
determine the waves transmitted and reflected for each incident wave. 

(The exact characteristic equations cannot be used since the derivatives 
become infinite at the discontinuity. ) Of course, an assumption regarding 
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the amount of energy added per unit mass flow must be made for the 
energy equation. This assumption is a very difficult assumption to 
make since the waves change the Reynolds number, velocity, and air- fuel 
ratio at the combustion zone. 

Two very simple assumptions were used for the two illustrative com- 
putations; namely, (l) the ratio of total temperature across the flame 
front was constant, and (2) the energy added per unit mass was constant. 
For an initial Mach number of 0.200 ahead of the flame and 0.400 behind 
it, these assumptions prescribed a constant stagnation-temperature ratio 
of 3*0482 for the first case and an increment in stagnation temperature 
of 2.0482 times the initial, undisturbed, stagnation temperature for the 
second case. A downstream-moving expansion was then assumed to strike 
the flame and theoretical computations made for the wave pattern due to 
the reflections and re- re flections from the flame discontinuity and 
nozzle. Under the limitation of a choked exit nozzle the waves rapidly 
damped out. The general form of the wave diagram is shown by the solid 
lines of figure 23 and the corresponding pressure-time curve at x = x^ * 

(located ahead of the flame by an amount equal to one-half the distance 
from the flame to the nozzle) for three different initial expansion 
strengths is plotted in figure 24. 

Computations were then made for a discontinuity in area in place of 
the flame front with the requirements that: (l) undisturbed conditions 

ahead of the discontinuity be identical to those for the heat addition, 
and (2) the medium-strength expansion have the same first upstream 
reflection for both the area change and heat addition. By trial and 
error the area ratio corresponding to the constant-stagnation-temperature 
ratio was found to be 0.5129 and that for the constant increment in stag- 
nation temperature to be O.38OI. The length of the system after the 
discontinuity in area change was determined so that an infinitesimal 
disturbance would require the same amount of time to travel from the 
discontinuity to the exit nozzle and back to the discontinuity for both 
the heat addition and contraction cases. In other words, if e a denotes 

the distance from the flame front to the nozzle and Cp the distance 

from the area change to the nozzle (see fig. 22), 
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The remainder of the wave pattern until* damping was completed was then 
computed, under the same restrictions as imposed on the heat-addition 
solution. Good agreement is shown in figure 24 for the pairs of 
pres sure- time curves for the corresponding conditions of (l) total 
temperature ratio equal to 3.0482 and area ratio of 0.5129, or ( 2 ) total 
temperature increment of 2.0482. and area ratio of O. 38 OI. For each of 
the three initial expansions., all the curves except that for the con- 
stant AX 0 will damp to the same final pressure (which is, of course,- 

different for each expansion) and return to' the original Mach number 
of 0.200 ahead of the discontinuity. Thus, for either of the heat- 
addition assumptions made, a corresponding area change could be employed. 

When the amplitude of the buzz in an actual ram jet is small, the 
variation of Reynolds number, velocity, and air- fuel ratio is also small. 
In a ram jet, however, there is also the added complication of coupling 
between the fuel- injection nozzles and the combustion zone whereby a 
wave originating from the combustion zone travels upstream' with a 
velocity a - u to the fuel-injection system and changes the air- fuel 
ratio of the fluid present at the fuel nozzles. This fluid then travels 
downstream with a velocity u to the combustion zone where the new air- 
fuel mixture bums' to generate a new wave. This phenomenon cannot be 
duplicated by waves whose speed is strictly limited to the local sonic 
velocity superimposed on the local fluid velocity. If the reflections 
from the shock at the inlet, from the nozzle exit, and from the waves 
striking the flame front, however, have a stronger influence on the 
pressure pattern than the waves generated by the change in air-fuel 
ratio, there is the possibility that the combustion in a particular 
buzzing ram jet may follow approximately a law of energy release which 
can be simulated by an area change. In that case the actual ram jet 
(fig. 22 ) with a finite combustion zone might be v tested for its behavior 
during buzz by the cold- flow test of the configuration (fig. 22 ) in which 
the combustion chamber is replaced by a shorter section with a gradual 
area reduction. 


CONCLUDING REMARKS 


Experimental instantaneous pressure records of simulated ram- jet 
models without heat addition showed that these models buzzed in a manner 
governed internally by quasi-one-dimensional (plane-wave) unsteady- flow 
theory, and not in accord with an analysis of the problem by the laws of 
Helmholtz resonator (spherical- wave) theory. A theory applicable to the 
buzzing problem was obtained through modification of the one-dimensional, 
unsteady- flow theory by introducing specified assumptions which greatly 
facilitated the computational efforts without impairing the usefulness 
of the theory for this flow. Because of the difficulty in computing the 
three-dimensional, unsteady, viscous flow when the shock was in motion 
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ahead of the cowling and in computing when the shock was in motion far 
from equilibrium with the waves moving upstream from the nozzle, it was 
necessary to use a few experimental points- at the start of the cycle as : 
a basis on which to compute the remaining pressure points of the cycle. 
Theoretical computations of the pressure- time curves, which were made 
on this basis for one supersonic-diffuser configuration with two differ- 
ent subsonic-duct lengths, agreed closely with the experimental records. 
The entrance and exit mass-flow - time curves, as well as the wave dia- 
grams, were also computed for these configurations by the unsteady- flow 
theory. ' 

The wave pattern causing the buzz was found to be a basic cycle 
which consisted of a downstream expansion reflected at the exit as an 
upstream expansion and a downstream compression reflected as an upstream 
compression plus a variable number of compression-expansion couplets 
following the basic cycle. The number of compression-expansion couplets 
depended on the valve-exit area and decreased as this area decreased] 
therefore, the frequency generally increased as the throttle was closed. 

The simulation of combustion in a ram jet by an area contraction 
for buzzing tests was discussed as a possible extension of the results. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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APPENDIX A 

DERIVATION OF CHARACTERISTIC EQUATIONS 


The general form - of the quasi- one-dimensional continuity equation 
is expressed in equation (54) in which the area A is considered to be 
a function of both time and position. With this equation, cases of a 
boundary layer growing with time as well as continued throttling may be 
treated: 


^-(pA) + i-(puA) = 0 
dt dt 


(54) 


Or, in another form, 


S log- A 3 logp 
— — - — + u — + 


St 


Sx 


A d log p 8 log- p du 

— + •—£ — + u — + — 

dt dx dx 


= 0 (55) 


The following two identities: 


a 2 = 7 RT 


j S _ _d_(heat added) 
R RT 


(56) 

(57) 


are now introduced. The quantity a in equation ( 56 ) is the speed of 
an isentropic wave of infinitesimal strength, or sound wave, which may 
move through a fluid of varying entropy. The "heat added" of equa- 
tion ( 57 ) includes the heat due to convection, dissipation, and sources. 
When the energy equation (equation ( 58 )) and the equation of state 
(equation ( 59 )) are combined, 


6 (heat added) = C y dT - dp 

P 2 


(58) 


p = pRT 


( 59 ) 
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with equations (56) and (57), the differential equations between the 
properties of state and the nondimens ional entropy are determined: 


2C 

d | = ~ y ~ d logg a - d logg p 


(60) 


c 2Cp ■ • 

d I = IT d log e a - d log e P 


' (61) 


Although these relationships were derived for a particular elemental 
volume of fluid, they also are true for all the various parts of the 
fluid and may he used when going. from one part of the fluid to another. 
This procedure is justified because the entropy is a property of state 

so that S(x,t) = s[p(x,t), p(x,t), a(x,t)j if all the fluid was origi- 

nally in the same state. Consequently, the following partial derivatives 
may be obtained: 


8 logg p 

dx 


6 logg p _ 

at 


6 loge 

Sx 


St 


2C V '8 log e 


4 


R 

ax 

Sx 

2C y 

a logg a 

.it 

R 

at 

St 

•2C p 

a logg a 

a s 

R 

R 

ax 

ax 

2Cp 

a logg a 

a £ 

R 

R 

St 

~ at 


(62) 


( 63 ) 


( 64 ) 


( 65 ) 


The momentum equation 


— + u — + = 0 

at a x pax 


( 66 ) 
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becomes equation ( 67 ) upon insertion of equations (5 6) and (59) 


Su + u — + — — Z = 
St Sx y Sx 


(67) 


Substitution of equation (64) into equation ( 67 ) and simplifying yields: 


Su Su 2 C V da 

— + u — + -r- a — 

St Sx R Sx 




( 68 ) 


since y = =£. 

W 

Similarly, substitution of equations (62) and ( 63 ) into equa- 
tion (55) results in: 


a 


S logg A 
St 


+ au 


S logg A 2 C y 
Sx R St 




Sa 

Sx 


au 



+ a 


^ = 0 
dx 


( 69 ) 


Addition of equations ( 68 ) and ( 69 ) and then rearrangement yields 


8 /2C V \ . \ 8 / \ 

^e- — - a + u + (u + a)~-( — 1 a + u 1 = 
8t \ R / ox\ R / 


/d log^ A 8 log A' 

-al “ — + u — — ] + 


8t 


8x 


y — a(tl + u I!) + & 

7 Vat a x / r 


it 

at 


6 | 

(u + a)— -5. 

dx . 

( 70 ) 


Similarly, subtraction of equation ( 68 ) from equation ( 69 ) and rearrange- 
ment gives 
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Uir a - u ) + (u - a - u ) * 


f b log A d log A' 


bt 


bx 


— ^ +a; ^, +L i 2 # +u y) + f 


y + < a 


- a) : 


( 71 ) 


The following identities are used: 


(1) Convective derivative — = ^ + u 

v ' Dt at ax 


(72) 


(2) Total derivative along a characteristic of slope = u ± a, 

CTG 


5 d / 4 _ \ d 

_ = + (u + a)— 


St dt 


bx 


(73) 


( 3 ) p = — a_ + _y_ = _a_ + 

y — \ 9* ° -no o 


*°oo a 0 ro 


R ^*Q ®*Q 

U 0O U 00 


(4) Q = 


- 2 a u _ a u 


7-1 a °oo a °» R a °oo a o« 


Insertion of these terms in equations (70) and (71) gives the character- 
istic equations in their -final form: 


8 P = a D lo 6 e A + 7-1 a D R + 1 _a_ ^_R 
St " a 0a> Dt 7 a 0 „ Dt 7 a 0 _ st 


SQ 

St 


OH _ a 


lo^A y i 


D I 


Dt 


R , 1 a 


8- 1 


7 a Q Dt 7 a 0 St 
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APPENDIX B 

COMPUTATION FOR CONFIGURATION I, = 44°, SUPERSONIC DIFFUSER 
WITH OVER-ALL RAM-JET LENGTH OF 1 % INCHES 


The first step in the computation for configuration I, 6 ^ = 44°, 

supersonic diffuser with over-all ram-jet length of 15 6 inches is the 
solution of the following equation (which is the same as equation (43) 
of the text) 




7+1 

2(7-1) 


For the effective area A e , the value of the area at x = 12.75 inches 
is chosen and. is equal to the area of the annulus formed by the inside of 
the pipe, which forms the principal part of the diffuser and by the sup- 
porting strut for the central body which extends to x at 15.5 inches 

(fig. 5) 


A e = |(2;062 2 -0.436 2 ) = 3 .I 9 O 


(72) 


Also, since to the order of accuracy of the "no- flow" shadowgraph 
measurements, -from figure 25 of- reference 3, the buzz should startat a 

V 

m„ 


mass-flow ratio 


= 1.000 because the conical shock falls on the 


cowling lip. When these , values are used with 7. = , 1.400 plus: the meas- 
ured- values of A^ = £.(]_, 2)^ square' inches, p^" = '87. 0 -pounds per square 
inch absolute (fig. 9)> P Q .=..101 .pounds per .square inch, absolute 

u oo 

(recorded during test), and the nozzle Mach number M^ = 1.90, one solves 


equation (43) for M e = 0. 1528; The effective entropy -^-(based on = 0 

is determined from the assunqrtion of steady" flow just prior to the start of 
buzz; thus, since the above data yield an effective stagnation pressure in 
the duct at start of buzz of 88.4 pounds per square inch absolute 


S„ 




Si 


Si : . • • 


^b _ ^e 
R R 


log’ ^^ : 0.1332J 


- ' : ( 73) 
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Similarly, 



_1 

2 

+ X O.1528 2 ) = 0.99767 

= 0.1528 X 0.99767 = 0.15244 


which yield starting values of the parameters P and Q 


P b = 5.1408 
Q b = 4.8359 


(74) 

(75) 


(76) 

(77) 


The reflection curve (solid line in fig. 25) for the entrance is 
computed as set forth in the theoretical section by adjusting the quasi- 
steady values of P 2 and by a factor which is a function of the 


difference 


£2 

R 


0.1332. Reflections at the nozzle exit are specified 


5 - M, 


by the relationship Qo = Po = 0.94069Pq (dashed line in fig. 25 ). 

-J 5 + o J 


The final curve (fig. 26) needed for the completion of the wave diagram 

_ \7 

is the pressure - sonic-speed curve Q -- — = 

87. 


te conq>j.eT.ion 03 

s_ . f s/a °- V 

r.o \o. 99767/ 


(see equation 21). 


Experimental values, designated by the flagged symbols on the pres- 
sure record (fig. 9) and wave diagram (fig. 8), are used to determine the 
downstream-moving waves when the normal shock is in rapid motion or ahead 
of the cowling. A linear interpolation of the variation of the param- 
eters P and Q is assumed for all waves. 


The remainder of the wave diagram is constructed by the two reflec- 
tion curves and then the pressures found on the third curve from the local 

sonic velocity ■ ~i~ 

a O 00 10 


Approximate values of the entropy at the entrance may be found from 
the quasi-steady flow to the inlet when the shock is ahead of the cowling. 
Consider the conditions existing (fig. 8) at t = 0.0193 a t the entrance 
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where P = 5.006l - and Q = 4.8359 so that =-0.0851, -2_ = 0.98420 

°oo a Ooo 

(in the imaginary lamina) and p = 79*10* The quasi- steady- flow energy 


equation gives the actual 


a _ 


= 0.999275 and, since at a free-stream 


Mach number of 1.90, p = 15*07 for p = 101.0, and 


= 0.76205, 


S 2 - Soo 


— — log^ °- 999275 _ x 19^0 = . 

R 7 - 1 ** 0.98420 e 15.07 


(78) 


The linearized approximation might he used instead to obtain the value 
S 2 

of IT 


s 2 - s «, 


So - S S - S 27 An 

= - - - + -2— — 22 = — + 0.1332 

R R R 7 - 1 a 


= 7 °‘9993 ~ 0*9842 + 0>1332 = 0.1058 + 0.1332 = 0.2390 (79) 

0*9993 


The ratio of the instantaneous entering mass flow to the mass flow at 
start of buzz is then, from equation (29): 


= /5.0061 + 4.8359x5.6061 - 4.8359 A _ 0*1058 \ = q>506 (8q) 

m b \5*l4o8 + 4.8359/5*1408 - 4.8359 \ 3.5 / 
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(a) Typical ram jet. 



(b) Hypothetical ram jet for one-dimensional, unsteady- flow 

analysis. 

Figure 1.- Sketch of typical ram jet and simplified one-dimensional 

counterpart. 
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(a) Variation of parameters across an entropy discontinuity 
where pressure and velocity are constant. 

Figure 3*- Diagrams showing how characteristic parameters P and Q 
may be referred to effective entropy by insertion of imaginary 
lamina of fluid in the flow field. 
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(b) Insertion of imaginary lamina of fluid (c) with entropy S^/R into 

flow (d) of different entropy returns parameters to valuo existing 
originally in state (3) to permit calculation of pressure and 
velocity.- * . 

Figure 3*- Concluded. 
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(a) Short diffuser, L = 34 inches. 



I i I I I l i 


(b) Medium diffuser, L = 156 inches. 



(c) Long diffuser, L = 390 inches. 


Gate valve 

(n) Location of pressure pickup, n (not necessarily same for various 
diffusers) 

r Location of 90° elbows 

0 .1 .2 .3 .4 .5 .6 .7 .8 

Position olong ramjet axis, x/L 

Figure 6.- Schematic diagram of ram- jet models tested showing: 

(1) relative length of central body to over-all length, and 

(2) locations of pressure pick-ups and .elbows. 
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Position along ramjet axis, x/L 


Figure 8.- Wave diagram for configuration I, 9 ^ = 44°, L = 156 inches. 

Flagged symbols denote points computed from experimental pressure 
record. 
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L- 72715 

Figure 10.- Low-speed, shadowgraphs during buzz of configuration I , 

= 44°, L = 156 inches. 
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Time , sec 


(a) Curves for gage 1 at - 0.074. 

Li 

Figure 11.- Theoretical pressure-time curves showing effect of variation 
m b 

Flagged symbols denote common pressures taken from 
experimental record. 
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Figure 16.- High-speed shadowgraphs during buzz of configuration I, 

0j = 44°, L = 34 inches. 
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Figure 18.- Pressure record for configuration I, Oj. = 44°, 
v L = 156 inches as exit throttle is rapidly opened. 
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Figure 18.- Concluded 
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Figure 20.- Pressure record for configuration II, 0^ = 43° 10’, 
L = 34 inches, showing unstable oscillation as buzzing begins. 



NACA RM L52A1 8 


CONFIDENTIAL 








KZj 1 

r so 

ECZ 1 

r-3, 

lv_- 

c 1 

C 132 

jr “ 03 

If— 

CZ 

ll— 

O ' 

jjjj^ 

0* 

d 06 

(CZ 

CZ 107 

ir" 

P 137 

CZ™ 

iT“ 

t 138 

mw 

lv__ 

ff “ 39 

«KmK 99 P^ 

cz ^ 

fcZ : 40 



L-72717 

Figure 21.- High-speed shadowgraphs during buzz of configuration II, 

Q z = 43° 10', L = 34 inches. 
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(a) Discontinuous heat addition in semi-infinite pipe. 



H — —M 

(b) Discontinuous area change in semi -infinite pipe. 



(c) Ram jet with finite combustion-chamber length. 



(d) Ram jet with gradual contraction of rear section. 

Figure 22.- Configuration pairs for possible . simulation of heat 

additions. 
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Distance along pipe f x 



Characteristics for heat addition 

— Characteristics for area change 

Figure 23.- Simplified wave diagram for heat addition and corresponding 
area contraction in semi-infinite pipe. 
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a 0(D t 

Time parameter , ^ 

Figure 2k.- Comparison of pressure-time curves for discontinuous heat 
addition and area contraction in semi- infinite pipe. 
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5.00 5.04 5.08 5.12 5.16 5.20 

Downstream characteristic parameter, or 

Figure 25.- Curves for determination of parameters reflected from 
cowling inlet and nozzle for configuration l f 9 ^ = 14°, 

L = 156 inches. 
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